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Abstract 

We study the two-dimensional contact process (CP) with quenched disorder 
(DCP), and determine the static critical exponents f3 and i'±. The dynamic 
behavior is incompatible with scaling, as applied to models (such as the pure 
CP) that have a continuous phase transition to an absorbing state. We find 
that the survival probability (starting with all sites occupied), for a finite- 
size system at critical, decays according to a power law, as does the off- 
critical density autocorrelation function. Thus the critical exponent z^||, which 
governs the relaxation time, is undefined, since the characteristic relaxation 
time is itself undefined. The logarithmic time-dependence found in recent 
simulations of the critical DCP [Moreira and Dickman, Phys. Rev. E54, 
R3090 (1996)] is further evidence of violation of scaling. A simple argument 
based on percolation cluster statistics yields a similar logarithmic evolution. 
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I. INTRODUCTION 



Phase transitions between an absorbing state (one admitting no further evolution), and 
an active regime occur in models of autocatalytic chemical reactions, epidemics, and trans- 
port in disordered media Paradigms of this sort of transition are the contact process 
(CP) and its simultaneous-update counterpart, directed percolation (DP) [^]. Since 
many-particle systems often incorporate frozen-in randomness, it is natural to investigate 
the effect of quenched disorder on an absorbing-state transition. Thus, some years ago, 
Noest observed that the critical behavior of disordered directed percolation is quite differ- 
ent from that of pure DP 0]. We recently studied the CP with quenched disorder in the 
form of random site dilution, and found logarithmic time-dependence at the critical point 
IP]. For example, P{t), the probability of survival, starting from a single active site, follows 
P ~ l/(lnt)" for large t. Such a form is incompatible with the scaling hypothesis that ap- 
plies quite generally to absorbing-state transitions 0]. This violation of scaling is consistent 
with Janssen's recent field-theoretic analysis, which shows that the resulting renormalization 
group equations have only runaway solutions Here we present further results bearing on 
the violation of scaling. In particular, we find that the exponent z/|| does not exist for the 
diluted contact process (DCP). We also propose an explanation for logarithmic behavior at 
the critical point. 

In the CP on the square lattice, each site is either vacant or occupied by a particle. 
Particles are created at vacant sites at rate An/4, where n is the number of occupied nearest- 
neighbors, and are annihilated at unit rate, independent of the surrounding configuration. 
The order parameter is the particle density p; it vanishes in the vacuum state, which is 
absorbing. As A is increased beyond Ac = 1.6488(1), there is a continuous phase transition 
from the vacuum to an active steady state; for A = A — Ac > 0, the stationary density 
p ~ A^. In the vicinity of the critical point the characteristic relaxation time r ~ lAI"*^!!^ 
and the correlation length diverges as ^ ~ |A|~'^^. In the DCP, a fraction x of the sites 
are diluted at random, and the birth-and-death process defining the CP is restricted to 
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non-diluted sites. (Further details on the DCP may be found in Ref. p.) 

To provide the necessary background we summarize the scaling behavior of the CP and 
allied models 0. Consider first the evolution from an initial configuration with just a single 
particle at the origin. The conditional probability of finding a particle at r, given that at 
time zero there was a particle at the origin, and that all other sites were vacant, obeys 



p(r, t; 0, 0) ~ ^''-'^"/2^(rV^^ At^/'^ii ). 



(1) 



Similarly the survival probability is expected to follow 



p{t) ~r^<i>(Ati/''ii). 



(2) 



{F and $ are scaling functions.) At critical (A = 0), equation (H) implies P(t) ~ t~^, while 
(H), when integrated over space, yields a mean population n(t) ~ t^. If we take the second 
moment (in space) of p{r, t; 0, 0) with A = 0, we obtain 



j.-q—dz/2 

R\t) ~ I x^f{x^/e) (fx oc e 



(3) 



For A > 0, the survival probability attains a finite asymptotic value: lim(_,oo P{t) = Poo ~ 
A^' 1^,1^. Several scaling relations can be derived, in particular, 

z = 2z/x/z/ii, (4) 

and 



(5) 



as well as the hyperscaling relation 



4:5 + 2r] = dz. 



(6) 



Eqs. (1) - (3) describe spreading from a single seed. Consider, on the other hand, a 
system at critical with all L'^ sites initially occupied, and let Pm{t) be the survival probability 
starting from this maximally occupied state. At short times Pm = 1 and the density is 



governed by the power law p{t) ~ t~^. Following this initial phase, Pm decays exponentially, 
with a characteristic time r ~ L^w/^^^ and the density in the surviving sample attains a 
quasi- stationary value ps ~ L~^^'^^. The scaling results summarized above have been amply 
confirmed for the CP and other models with a unique absorbing configuration, and have 
been extended to models possessing multiple absorbing configurations P,pi|. 

In our recent study of the DCP, we found that P(t), n(t) and -R^(t) display logarithmic 
time dependence at critical, which is incompatible with the scaling forms, Eqs. (1) and (2), 
describing the pure model. In this work we report further results on static and dynamic 
critical properties, in an effort to determine the extent of the scaling violation, and to 
understand its origin. In Sec. II we analyze the stationary density (in the supercritical 
regime), and the quasi- stationary density (at critical), to obtain estimates of f3 and z/_|_, and 
also examine the stationary density- density correlation function, which decays algebraically. 
Sec. Ill concerns the survival probability, Pm, starting from a maximally occupied state; it 
also decays algebraically, roughly as predicted by a simple probabilistic picture. Another 
simple argument is presented in Sec. IV, for the logarithmic decay of the survival probability. 
We conclude, in Sec. V, with a discussion of our main results, and of the reason for violation 
of dynamic scaling in the DCP. 

II. STATIC BEHAVIOR 

A. Stationary Density 

We determined the stationary density p for dilutions x = 0.05, 0.1, 0.2, 0.3 and 0.35, by 
the following simulation procedure. After generating a disorder configuration on a lattice of 
L X L sites, and initially occupying all non-diluted sites, we permit the system to relax for a 
time tR, and then accumulate data on the density for a period of duration ts- This process 
is repeated Nt times (with a new disorder configuration for each trial), and the density 
pit), computed over the Ns trials that survive up to time + ts, is examined to verify 
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that sufficient time has been allowed for relaxation. (If not, is increased accordingly.) 
We then take the mean density for each trial during the observation time < t < ts, and 
compute the mean and standard deviation over the sample of Ns independent trials. We 
used lattice sizes L = 100, 200, 400, and 800, increasing L as we approached the critical 
point, and checking (for L < 400) that estimates for the stationary density agreed for two 
lattice sizes. The sample size Ns ranged from 50 to 100. We employed tn = ts = 5 x 10^ - 
2.5 X 10^, the larger values reflecting slower relaxation near Xc{x). (We use the critical point 
estimates \c{x) obtained via time-dependent simulations in Ref. p.) 

Our results for the stationary density, shown in Fig. 1, reveal a crossover between the 
DP value of /5 ~ 0.58 at small dilutions (x = 0.05 and 0.1), and a larger exponent as one 
approaches Ac. The data for x > 0.1 yield exponents in the range 0.89 - 0.98, suggesting 
that P = 0.93(5). (Figures in parentheses denote statistical uncertainties — one standard 
deviation.) Similarly, an analysis of the stationary mean-square density supports ~ 
with 2/5 = 1.84(4). Thus the order parameter exponent agrees, to within uncertainty, with 
our earlier result /?' = 0.99(3) for the exponent governing the ultimate survival probability 
1^. (For the pure CP, /5' = /5 ^ 0.58.) While our estimate is not far from Noest's result, 
P = 1.10(5), we regard it as excluding that value; we do not feel, however, that our data 
rule out P = 1. 

B. Quasi-Stationary Density 

Rather than studying spatial correlations directly, we determine i>± by analyzing the 
quasi- stationary density ps at critical as a function of the system size; this yields P/i^±, as 
noted in Sec. I. We studied lattice sizes varying from L = 8 to L = 128, averaging over 
2 X 10^ to 10^ independent runs, of duration tg = 10^ to tg = 10^. (The larger sample sizes 
and longer run times apply to the larger L values.) We show in Fig. 2 a log- log plot of 
Ps{Xc,L) versus L for dilutions ranging from 0.02 to 0.35, along with the slopes of linear 
least-squares fits to the data for L > 16. (The uncertainty in the slope ranges from 1 to 
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3%.) For X < 0.1, the scaling of ps is similar to that found in the pure CP, for which 
/3/z/_L = 0.80(3). But for larger dilutions we see a steady increase in the slope; we estimate 
/3/z/_L = 0.93(3) for the DCP. Combining this with our estimate (3 = 0.93(5), we obtain 
= 1.00(9) for the DCP. This is significantly lower than, but still consistent with, Noest's 
earlier estimate of 1.17(10) On the other hand, a theorem of Chayes et al. requires 



di'± > 2 for models with quenched disorder, or u± > 1 here [|r^; the upper range of our 
estimate is consistent with this result. (Note that if we use our earlier estimate, /?' = 0.99(3) 
in place of our /3 value, we obtain z/^ = 1.06(7).) 



C. Stationary Density-Density Correlation Function 

The stationary density- density correlation function is defined by 

C{t) = {a,{to + t)aM) , (7) 

where ai{t) = 1 (0) if site i is occupied (vacant) at time t, and the average is over realizations 
of the process and over disorder. It is understood that to is sufficiently large that the r.h.s. 
is independent of to. For the undiluted CP, the r.h.s. is then independent of i as well, and 
AC = C{t) —J)'^ exp(— t/r) for large t. r is a characteristic relaxation time diverging as 
r ~ A^^'ii near the critical point. 

We studied the density-density correlation function in the DCP at x = 0.1 and 0.3, 
evaluating the r.h.s. of Eq. (|^) for a single site i (the first non-diluted site to be generated), 
on L X L lattices with periodic boundaries, using 500 — 2000 independent realizations of 
the disorder. Fig. 3 (inset) shows a typical evolution, with C{t) slowly approaching an 
asymptote. Coo- The main graph shows that the excess, AC = C{t) — Coo ~ t~^, so it 
cannot be characterized by a relaxation time. The exponent b varies from about 0.7 - 0.8 
well above Ac, to about 0.5 near Ac. Note that for & < 1, even the alternative expression 

roo 

T = / dtACit) , (8) 
Jo 
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similar to that employed by Noest, is undefined. We also observe a power-law approach 
(with an exponent of about 1/2) of the critical quasi-stationary density p^, discussed in the 
preceding subsection, to its asymptotic value. 

III. DYNAMIC BEHAVIOR 

Consider the pure CP on a lattice of L"^ sites, starting with all sites occupied. Since we 
are dealing with a finite system, there is a well-defined lifetime t( A, L) and at long times 
the survival probability Pm{t) exp[— i/T(A, L)] for large t. Just at the critical point the 
lifetime has a power-law dependence on L: r{0,L) ~ L'^ii/'^-l. For A > 0, vacating {L/^Y 
independent regions simultaneously is an exponentially rare event, and we expect the scaling 
form 

t(A, L) - L"ii/"^ exp[c(LA'^^)'^], (9) 

(c is a constant), as is confirmed by the data shown in Fig. 4. (The data also appear to 
scale for A < 0, but with a different scaling function.) 

We studied the relaxation from a maximally occupied state (all non-diluted sites occu- 
pied) in the critical DCP. Figs. 5 and 6 show Pm{t) for x = 0.1 and 0.3, respectively, for 
L = 32, 64, and 128. (We studied samples of 5000, 2000, and 1000 trials for L = 32, 64, 
and 128, resp.) From the figures it appears that following the initial stage, the survival 
probability crosses over to a nonuniversal power law, with an exponent that decreases with 
L. (For X = 0.1 we find ~ t"" with a = 1.8, 1.1, and 0.86 for L = 32, 64, and 128, resp.; 
for X — 0.3 the corresponding powers are 0.71, 0.55, and 0.37.) Since the asymptotic decay 
of Pm at critical is nonexponential, there is no characteristic lifetime for the process. The 
initial stage, during which Pm = 1, is characterized by a correlation length ^(t) < L. During 
this phase we find rim ~ with 5 ~ 0.47, as for the pure model. Once Pm starts to decay, 
Um crosses over to a different, nonuniversal power law, as seen in Fig. 7. 

We note in passing that while the Pm data are certainly inconsistent with exponential de- 
cay, a slight (downward) curvature suggests a faster than power- law decay. In fact, somewhat 
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more linear plots are obtained using the form InP^ oc (Int)"'^, with a in the range 1.5 — 3 
depending on the data set. The data are too noisy to permit a definite conclusion regarding 
the form of Pm-, but it is interesting that a simple argument yields the modified power law. 
We begin by supposing that because of fiuctuations in the disorder, the probability of a 
lattice of L'^ sites having an effective distance from criticality of A is -P(A) ~ exp(— foL'^A^). 
(Here we ignore the spatial inhomogeneity of the disorder, and invoke the central limit the- 
orem.) If we then assume that such a lattice is characterized by a lifetime r( A, L) scaling 
as in the pure model, Eq. (||), the survival probabihty is given by 



Maximizing the argument of the exponential to extract the leading behavior at large t, we 
obtain InP ~ {Intfl'^''^. (The exponent ~ 1.37 for d = 2.) 

This argument, which treats fluctuations in the disorder as if they were spatially homo- 
geneous, and uses pure-model scaling, is clearly inadequate to deal with the true subtlety of 
the DCP. That it may yet contain some germ of truth is suggested by our flnding that when 
we do not average over disorder, the decay of Pm is exponential, as in the pure CP, but with 
a lifetime particular to the disorder set generated. (One is naturally interested in knowing 
the distribution of the relaxation time. This poses a formidable numerical task that we hope 
to address in future work.) In summary, the relaxation of the DCP from a maximally occu- 
pied state is similar to that of the pure model during the initial stage, in which correlations 
have yet to grow to the size of the system. But afterward the evolution follows nonuniversal 
power laws (or modifled power laws) and we cannot deflne a characteristic lifetime. 



In this section we propose a simple explanation for how logarithmic time-dependence 
arises in the critical DCP. Consider the survival probability P{t) starting from a single 
occupied site or seed located at the origin, O. Clearly, the trials that contribute to P{t) 




(10) 



IV. CRITICAL DYNAMICS 
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at large t are those in which the seed happens to fall in a large, favorable region. To 
make the notion of a "favorable region" somewhat more precise, imagine taking the disorder 
configuration on a cube of L''- sites, filling space with periodic copies, and running the contact 
process with A = Ac(x) on this lattice. For some disorder configurations — the favorable 
ones — the process will in fact be supercritical, because the fraction of diluted sites is < x, 
or because of a particularly advantageous arrangement of the diluted sites. Such regions are 
characterized by an "effective distance from criticahty" Ag// > 0. Any disorder configuration 
contains both favorable and unfavorable regions. If the seed lies in a favorable region, we 
can define its "domain" as the maximal connected favorable region containing O. (The 
domain is surrounded by unfavorable regions, which impede the spread of the process.) For 
simplicity, we suppose that on a domain of V sites, having some Ag// > 0, the DCP has a 
lifetime r ~ explcVAfj^j], as it would on a compact region, in the pure model. From the 
central limit theorem, the typical value of Ag// on a domain of V sites ~ F"^/^, yielding a 
lifetime r ~ exp[cV^~''''^-^/'^]. At time t, only trials whose seeds happen to fall in a domain 
with T >t, or V > c'(lnt)2/(2"''''^) survive. 

It remains to estimate the probability p{V) that O belongs to a domain of V sites. To 
do this, note that a domain is a kind of percolation cluster. The precise definition of the 
sites in this percolation problem is unclear (we might imagine averaging over small blocks 
of sites in the original lattice), as is the connectivity rule (next-nearest neighbor blocks, for 
example, might effectively be connected). But we should expect the associated percolation 
model to be isotropic and of finite range. Moreover, domain percolation must be critical at 
Xc{x). If it were supercritical, the contact process would be able to spread into an unbounded 
domain, and so would itself be supercritical. Similarly, if domains were subcritical, their 
size distribution would decay exponentially, and the CP would be subcritical. At critical, 
the domain size is power-law distributed: p{V) ~ for large V, with Tp the usual 

percolation cluster-size exponent {Tp — 187/91 in two dimensions). Combining this result 
with the lifetime estimate, we have 
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P{t) ~ / ~ /l^^)-2(r,-2)/(2-d.,)_ 

As in the argument (Sec. Ill) for the survival probabihty starting from a maximally occupied 
state, the effect of inhomogeneity in the disorder is greatly oversimplified. Inserting the 
known values of and of Tp, we obtain P{t) ~ (lnt)~°-^, whereas the exponent we observed 
in simulations ||^ is much larger, and nonuniversal, ranging from about 2.7 at a; = 0.35 to 
4.6 at X = 0.1. Thus we offer the above argument without any claim of quantitative validity, 
but rather to show how a simple treatment of disorder leads naturally to logarithmic time- 
dependence, and in the hope that it may form the basis for a more convincing approach. 

The probabilistic arguments suggest that it may be possible to understand how anoma- 
lous dynamics arises from an average over disorder. Here it is important to recall Noest's 
analysis of contributions from exponentially rare, favorable disorder configurations to the 
survival probability P{t). By deriving upper and lower bounds on the survival probability, 
he was able to prove power-law decay of P{t) in a Griffiths phase for A > Ac(0), in the non- 



percolating regime [|T2[- Arguments of a somewhat similar nature were advanced by Bray in 
his discussion of the relaxation of diluted spin models ||13 . 



V. DISCUSSION 

We have found that some aspects of the diluted contact process exhibit the same sort 
of critical behavior — albeit with different exponents — as seen in the pure model. Other 
features — spreading from a seed, the density-density autocorrelation function, C(t), and 
the survival probability Pm{t), starting from a maximally occupied state — do not follow 
the usual scaling, and are nonuniversal. The anomalous properties are all connected with 
dynamics: normally, C(t) and Prn{t) decay exponentially, with the diverging lifetime serving 
to define the exponent v\\ through r ~ A^'^n. Here no such definition is possible. Consistent 
with this, the spreading exponents 5, and r/, which are formally zero, are connected to 
v\\ via the scaling relations Eqs. (4) — (6). Since (3 and vx. are in fact finite, Eqs. (4) 
and (5) suggest that u\\ is infinite. Given Janssen's recent results it is of interest to 
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know whether our finding of power-law static behavior, but anomalous time-dependence, is 
compatible with a field-theoretical analysis . 

Some insight into the violation of dynamic scaling may be gained by returning to Eq. 
(y): the exponential dependence of the lifetime (in the pure CP) upon L and A suggests 
an extreme sensitivity of dynamic behavior to disorder. We might expect dynamics to be 
dominated by the extremes of local fiuctuations in the disorder. The result is that a property 
such as the relaxation time r for the density- density autocorrelation function is non-self- 



averaging, i.e., it does not converge to a limiting value even as L ^ oo |14|J15|| . Another 
manifestation of non-self-averaging is the dependence of P{t) and even at long times, 
on the location of the seed at time zero. These features are dominated by local fiuctuations, 
rather than by the properties of a "typical" disorder configuration. This in turn suggests that 
further insight may be gained by studying critical behavior for fixed disorder configurations, 
in order to determine the statistical distributions of various system properties, and of the 
domains defined in Sec. IV. In this way the primitive probabilistic arguments for and 
P{t) could be honed into a quantitative description of how anomalous behavior emerges in 
the average over disorder. 
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Figure Captions 

FIG. 1. Stationary density, p, versus A = A — Ac in the diluted contact process, x: x — 0.05; 
filled squares: x — 0.10; o: x — 0.20; O: x — 0.30; •: x — 0.35. Figures denote the slopes of 
the various straight lines. 

FIG. 2. Quasi- stationary density ps{Xc,L) versus L, for dilutions x = 0.02 (filled squares); 
0.05 (A); 0.10 (x); 0.20 (o); 0.30 (o) and 0.35 (□). The lines are linear fits to the last four 
data points in each set; figures indicate the slope. 

FIG. 3. Main graph: Excess density-density correlation function versus time in the DCP for 
X = 0.3 and A = 2.60 (top), 2.70 (middle), and 2.47 (bottom). The inset shows C versus t 
for X = 0.3 and A = 2.70. 

FIG. 4. Semilogarithmic scaling plot of the lifetime in the pure CP. Squares: |A| = 0.01; 
diamonds: |A| = 0.02; •: |A| = 0.05. 

FIG. 5. Survival probability Pm versus time in the critical DCP starting from a maximally 
occupied state, for dilution x = 0.1 and A = Ac = 1.8464. +: L = 32; o: L = 64; o: L = 128. 

FIG. 6. Same as Fig. 5, but for x = 0.3, A = Ac = 2.47. 

FIG. 7. Mean number of particles in the critical DCP with x — 0.3. Symbols as in Fig. 
5. 
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